Dynamical Structure Factor and Spin-Density Separation for a Weakly-Interacting 

Two-Component Bose Gas 
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We show that spin-density separation in a Bose gas is not restricted to ID but also occurs in 
higher dimension. The ratio (a) of the intra-species atom-atom interaction strength to the inter- 
species interaction strength, strongly influences the dynamics of spin-density separation and the 
elementary excitations. The density wave is phonon-like for all values of a. For a < 1, spin wave is 
also phonon-like. The spin waves have a quadratic dispersion in the a = 1 coupling regime, while 
in the phase separated regime (a > 1) the spin waves are found to be damped. The dynamical 
structure factor (DSF) reveals two distinct peaks corresponding to the density and spin waves for 
a < 1. For a > 1 there is only one DSF peak corresponding to the density wave. 
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Spin-density (charge) separation is a remarkable fea- 
ture predicted for one-dimensional interacting spin-1/2 
fermions [l[ and widely searched in condensed matter 
systems. Consequently, its investigation in atomic sys- 
tems could be of interest to different fields of physics. 
Unlike the higher-dimensional fermionic systems, where 
elementary excitations normally carries both spin and 
density(charge) degree of freedom, the collective excita- 
tions of the one-dimensional Fermi system separate into 
two distinct modes, spin and density waves due to the 
fact that the interaction in one-dimensional system lead 
to a Luttinger liquid state with bosonic excitations 
This behaviour is a hallmark of collective effects caused 
by interactions. For bosons, using two-components, cor- 
responding to two hyperfine states of cold atoms [2| al- 
lows us to study the (iso)spin waves as the relative spa- 
tial oscillations of the two-components. Till now, the 
study of spin-density separation has been limited to one- 
dimensional fermions Q and bosons 0, @ since, both 
these systems are believed to belong to the same Lut- 
tinger liquid class, which leads to the spin-density sepa- 
ration. Contrary to expectations, in this Letter we will 
show that for bosons spin-density separation not only 
exists in one-dimensional systems but also in higher di- 
mensions since for bosons, excitations are always collec- 
tive excitations in all dimensions due to the presence of 
the (quasi)condensate fraction and that the Luttinger liq- 
uid approach is not essential to describe the spin-density 
separation in bosons. We also compute the dynamical 
structure factor which reveals distinct features of spin 
density separation in all dimensions. 

We start with the Lagrangian density for two- 
component Bose gas at zero temperature: 
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real time, and here we set h = 1. Also \ii , m = \<Pi\ is 
the chemical potential and the particle density of the i th 
component and gd^j > is the repulsive effective atom- 
atom interaction between the i th and j th components. 
Here we consider bosonic atoms of the same isotope of 
mass m but having different internal spin states, there- 
fore we have g d ,n = gd.ii = 9d and g d . 12 = g d ,2i = 9d- 
For simplicity, we consider the same average atom num- 
ber density for the two components, i.e. fi\ = n 2 = h. 
For a 3D Bose gas , gs = ^ita^/ m d, 0). Here 03 is 
the 3D scattering length. For lower dimensional Bose 
gas in a 3D trap with longitudinal harmonic trapping 
frequency uj_, g 2 = An /(m In na 2 ) [1, @, E3] with the 2D 

scattering length given as a 2 ~ 7Ale~y /7r/a ^ muj± [n| and 
gi = 2oj±a,3 with 03 <C l/^/muj± [lij ]. The behavior of 
the system depends crucially on the dimensionless pa- 
rameter 7d(7d) = mgdn 1 ~ 2 l d {mg d n 1 ~ 2 / d ). For the gas 
to be weakly interacting, we must have 7^(7^) <C 1. The 
chemical potential fii of the i th component is determined 
by the condition J2ij 9d.i 3 nj = /^n 4 . 

To understand the low-energy excitations in two- 
component Bose gas, one can derive a low-energy ef- 
fective hydrodynamical Lagrangian that contains only 
modes related to the low-energy excitations (g, LL3J • We 
write the Boson field ifi in the terms of the number den- 
sity rii and the phase Oi as ifi = nie l6i . In the weak cou- 
pling regime the phase changes slowly in space while the 
density fluctuates fastQ, therefore one can integrate out 
the high energy fast density fluctuation [3] to obtain the 
effective hydrodynamic action. We introduce the density 
fluctuation drii as rii = n + 5rii. In terms of the new ba- 
sis, 5n p ( a) = (Sni ± 6n 2 )/V2 and 6 p{a) = (6>i ± 6> 2 )/V2, 
the action obtained from the Lagrangian density can be 
rewritten as 



where tfi — (fi(r,t),i = 1,2 is the field representing two 
different Bose particles, r is the space coordinate, t is the 
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where n p (n <7 ) = v2n(0) and 
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FIG. 1: Dynamic structure factor for all dimensions in the 
hydrodynamic regime and at zero temperature for a = 0.5 
and q = 1.25, shows two distinct peaks corresponding to the 
density and the spin waves, centered at v p q and v a q, respec- 
tively. The one dimensional structure factor is found to be 
a delta function, while the two- and three dimensional DSF 
is broadened because of Beliaev damping. In two dimension 
the DSF for the density wave is broader compared to that of 
spin waves, while the width of the three-dimensional peaks 
remains the same. There is no specific reason why the width 
of the 2D peaks is larger than both the ID and 3D peaks, 
this depending on the choice of the parameters . 
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For a < 1, after performing two Gaussian integrals, 
the effective action has the form[13] 

5 cff = f d d xdt ]T ^ (|9 4 0a| 2 - ^|V r A | 2 ) , (3) 
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where <^ (<7) = e i6 «M, X P (<r)(= l /9d,p (a)) is the density 
(spin) compressibility and v p / cr \(= \fngd,p(<r)/ m ) is the 
sound velocity of the density(spin) mode. Here we as- 
sumed that the fields dp(a-) vary slowly in space and we 
v 2 

have dropped the g^r- term. The effective action ([3)) 
describes the low-energy excitations of two sound waves 
with linear dispersions u> p ( a ) = v p t a \k. The bosons split 
into two gapless modes, namely density mode and spin 
mode, propagating with different velocities. The den- 
sity wave propagates faster than the spin wave, which 
can be seen by the relation v a /v p » — a)/ (I + a). 
In this regime the energy gap of the lowest excitation 
above the ground state is zero. Such systems have a di- 
verging length scale determing the exponential decay of 
equal time correlatons in the ground state, which defines 
the quantum critical behavior. Therefore the systems for 
a < 1 lie at the quantum critical points 

The meaning of the low-energy effective Lagrangian 
([3]) is that the bosonic system separates into two inde- 
pendent degree of freedom, i.e. spin and density. Unlike 
in fermionic one-dimensional systems, we do not need 
the bosonization method to obtain the spin-density sep- 
aration, the only thing we need is the (quasi) condensate 
density n to have fluctuations around it. This can 



be fulfilled in all dimensions at zero temperature for 
bosonic systems. The Bogoliubov energy dispersion re- 
lation of one-component interacting Bose gas is e(k) = 
y/ ((fc 2 /2m) 2 + g d nk 2 /m) [lij]. For the two component 
Bose gas, replacing the interaction g d with gd, P (a)i we 
obtain two branches of the excitations 



W*) - V( fc2 / 2m ) 2 + gw(1±a) " fc 2 , (4) 
V m 

which is in agreement with the result obtained by the 
semiclassical method (l6| . From the dispersion relations 
(HI we can define the chemical potential for the density 
and spin waves as /i p ( CT ) = gd.p(a)^- 

For a — l(gd — g'd)j only one Gaussian integral can 
be performed in action {2} giving the gapless density 
wave with linear dispersion. However, one obtains a 
quadratic dispersion for the spin-wave excitations, in 
agreement with 517(2) symmetry [j| [l^. This effect 
can also be seen from the Bogoliubov excitations e p — 
\J (fc 2 /2m) 2 + 2gnk 2 /m and e CT = k 2 /2m by replacing 
9d. P = 2gd and gd, a = in J4}. In this case, due to the 
SU(2) symmetry, the eigenstates are classified according 
to their total spin 5 ranging from to N/2, and according 
to recent result by Eisenberg and Lieb [18|, the ground 
state is fully polarized (5 = N/2). In one dimension, the 
ground state is described by Lieb-Liniger(LL) model of 
one-component interacting Bose gas [19j . for which the 
elementary excitations in the weak-coupling regime are 
density waves [20j j . and the system is ferromagnetic. 

In the case of a > l(gd < 3^), we found gd.a < 0. This 
implies, v„ (= Jj^i/m) in the long wave length limit 
is imaginary. The spin waves become unstable and damp 
out in the thermodynamic limit. Therefore we obtain a 
phase separation of the two-component Bose gas [3| • 

The dynamical structure factor (DSF) of many-body 
system is defined as follows 

S p[ *)(q,u>)= J d d xdte J ^ t ^ x )( ( 5n pW (x,<) ( 5n p(ff) (O,0)), 

(5) 

where (• • • ) can be calculated using path integral with 
the effective action. Experimentally, one can measure the 
dynamical structure factor using Bragg spectroscopy 24 1 . 

For a < 1, from the action [2l one can get 
the equation of motion for 5n p ^ as 8n p ( a ) (x, t) = 
~^/dd,p{a)dt0p(a)(' x -,t)- From the quadratic Lagrangian 
density, the DSF can be obtained as 



Sp(a)(<l,u) 



Im 



X P (a)V 2 p(a) q 2 



(?)' 



(6) 



where oj p ^(q) = v pi^q+iT p ua{(}) with the quasiparticle 
decay rate T p{a ){q). 

In order to obtain the DSF, one has to find the com- 
pressibility Xp(cr): velocity fp^) and decay rate T p ^(q) 
in terms of the dimensionless parameters 7<j )P ( cr ). Using 
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the macroscopic argument, the compressibility Xp(o-) is 
related to the energy E(g p(a)l n) as x~^) = V7f£ with 
the constant system size: V = L d and density: n = N/V. 
Similarly, the sound velocity can also be obtained using 
the macroscopic energy spectrum as v p ( a ) = (^fpf) 1 ^ 2 



with constant particle number N 21| . The way to obtain 
the ground state energy spectrum is diverse and depends 
on the dimension. As indicated by Beliaev 22j, the di- 
mensional dependent decay rate r p ( CT ) (q) is caused by the 
process of a long wave-length phonon decaying into two 
phonons and it can be calculated for small momenta us- 
ing the formula § 
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J d d fc|q||k||q-k| 
-e p (<T)( k ) -e p (o-)(q-k)). 



(7) 



In 3D, the dimensionless parameter 73 = Ana^n 1 / 3 . 
In this case, the requirement for a dilute gas na 3 -C 1 
corresponds to the weak-coupling condition 73 -C 1. The 
ground state energy was given for the first time by Lee 
et al. @ as E = Nn 2 / 3 /{2m) l3 (l + 16^ /2 /5ir 2 ). The 

ground state compressibility and velocity are given by 

3 
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It^^s p(<t)) 2 ' respectively. The decay rate for 3D system is 
obtained from eq®: T p{<j) {q) = T(q) = g^L- We 
can see that the decay rates for density and spin waves 
are equal and proportional to q 5 . The DSF for uj > can 
be approximated as 



Xp(a)Vp(a)q^p(a){q) 



(u-v p(a) q) 2 +Y 2 p{a) {q) 



(8) 



In the Bragg scattering experiment, one should obtain 
two peaks centered at v p i a \q for the cross section with 
the width T(q). 

For 2D Bose gas,renormalization-group analysis 
231 ] shows that the interaction of the 2D dilute gas is 
marginally irrelevant only in a dilute limit specified by 
In In 72 3> 1. The corresponding ground state energy for 
a weak-interacting gas is given by E = iVn/(2m)72(l — 
C72) where constant C <C 1 is not universal but model- 
dependent due to the marginal intcraction[23]. The com- 
pressibility and velocity for spin and density-wave excita- 
tions are = 52, P (<r)(l ~ {C - ^)l 2 , P (a)) and v p{a) = 

(1 - (C- ^r)72,p(<T)) 1/2 ■ The Belieav decay rate 



"72,p(<r) 



can be obtained by the integral (J7J): T p r a -\ (q) = 
Therefore the DSF §8$ has a broader width for density 
waves than spin waves. 

In the case of one dimension, contrary to 2D and 3D 
systems, the weak coupling means that the system is 
in the high density regime because 71 = mgi/n. In 
this regime, Lieb and Liniger [l9| first gave the ground 




FIG. 2: Dynamic structure factor for 5(7(2) in symmetric 
Hamiltonian in 3D. Note that uj is in the units of gdn and 
q is in the units of ^/mgdn. The DSF of the density waves 
varies linearly with q, while the DSF of the spin waves shows 
the quadratic dependence on q. 



state energy as : E = ^-71 (1 ~ SfVtT)- A few al g e " 



bra leads to the compressibility and sound velocities as 

Xp£) = -SfVTTJw) and v pW) = \f^¥^( 1 ~ 

^^/Ti.pfcr)) 1 ^ 2 - For ID, one obtains no decay rate. The 
reason is that the scenario for one phonon decaying into 
two phonons cannot exist due to the fact that energy con- 
servation law in eq® cannot be fulfilled in ID. Therefore 
two sharp peaks should be observed in the Bragg scat- 
tering experiments. Figure [T] illustrates and summarizes 
the results obtained above for the DSF in all the three 
dimensions. 

In the case of a = 1, the situation changes. For the 
density waves the dynamic structure factor remains the 
same as that in two-sound regime, while the DSF for 
spin-wave exitation alters due to the dramatic changing 
of the dispersion from linear to quadratic. In order to 
calculate the DSF one can use the effective Hamiltonian 
in the weak-coupling regime: 



2^ e P b l b p 




(9) 



: p 2 /2m, the Bo- 
e 2 + 2/idnep [H and the chem- 
ical potential : fi^ = 2g^n . Using 8n a — \/E(tf + b), the 
DSF can be related to the imaginary part of the Green 
function as S(q, uj) — n!mG(q, uj) where G(q, uj) is the 
singe particle Green function of the spin operators b q and 



with the spectrum of free spin waves e p 
goliubov spectrum e p - 
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6j. Therefore the DSF for the spin waves reads 



S(q,u) = 



(10) 



where the effective mass m* is determined by the equa- 
tion: m/m* = (1 + 2/md 2 ReY,(p)/dp 2 )(p = 0) with the 
self energy defined as S = G~ 1 (gd) — G~ 1 (g c i = 0), 
and the decay rate: Td, a (q) = ImE(g). To the sec- 
ond order diagram for the self energy S, one obtains 
the inverse effective mass related to the dimensionless 
parameter 7^ = ^n" 2 ^ as mjm* d = 1 — OLd^J 2 with 
ad = 2/37T, 1/2-7T, 1/8-7T for one, two and three dimen- 
sions, respectively. The decay process depends on the 
spin-phonon interaction which requires the energy con- 
servation: e g _fc + efe = e q with the spin momentum q 
and phonon momentum k. For q < ^Jm/Ad, this condi- 
tion cannot be fulfilled, therefor T da — 0, i.e. S(q,u>) = 
n5(ui — ujq). For q > y/mfid , an aproximation can be 
obtained as follows : V da (^/rfTJI2(l + 5)) = PdfJ-d^ 3 ^ 2 ^ 1 ^ 1 
for 8 < 1 with d = 0, 1/8tt, 2/3vr for d = 1,2, 3, respec- 
tively. The eqn. (flQ|) shows the fact that the excitations 
for spin waves for the Hamitonian with SU (2) symme- 
try are not sound-like, but particle-like, with the DSF 
centered at the position proportional to q 2 instead of q. 
Similar to the two-sound mode regime (a < 1), the DSF 
for one dimension is a delta function due to the fact that 
the energy conservation relation for a particle emitting 
a phonon cannot be fulfil in one dimension at zero tem- 
perature. Fig. [2] shows the linear dispersion of the DSF 
for the density waves and the quadratic dispersion of the 
DSF for the spin waves. The delta function behavior are 
shown for low momenta. 

For phase-separation regime (a > 1), the spin waves 
are thermodynamicall unstable, therefore only density 
waves exists. The dynamic structure factor of spin waves 
is smeared out and there exists only one peak in the DSF, 
which is different from the other regimes. This property 
can be a prominent signature for checking whether the 
system is phase separated or not. 

In this Letter we have shown that, unlike fcrmionic sys- 
tems, spin-density separation in two-component bosonic 
system is a more generic feature and occurs in all dimen- 
sions. The density wave is found to be phonon like for 
all dimensions and coupling regimes. However, the spin 
waves in the two-sound regime (a < 1) show a linear 
dispersion (phonon like) and the DSF for all dimensions 
show two distinct peaks corresponding to the density and 
the spin waves, centered at v p q and v a q, respectively. 
In the same regime, the one-dimensional structure fac- 
tors are found to be delta functions, while the two- and 
three-dimensional DSF is broadened because of Beliaev 
damping. The spin waves show a quadratic dispersion in 
the SU(2) symmetric regime (a — 1) and the DSF for 
all dimensions also show two distinct peaks centered at 
v p q and q 2 /2m* { . The spin wave is damped in the phase 



separated regime (a > 1) and there is only one peak cor- 
responding to the density wave. These are interesting 
signatures of spin-density separation to look for using 
Bragg spectroscopy where, the response of the conden- 
sate to a two-photon Bragg pulse is measured 24]. The 
difference between spin and charge velocities allows us 
to have spin and charge wavepackets moving at differ- 
ent velocities. An optical potential generated by a laser 
tuned, e.g., between fine-structure levels of excited alkali 
states transfers momentum solely to the the spin waves, 
while an optical potential far detuned will act solely on 
the density waves Q. One can also coherently excite 
the spin waves and the density waves simultaneously and 
then probe the two waves with a second laser pulse at a 
later time. Spin-charge separation manifests itself in a 
spatial separation of the spin and density wavepackets. 
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